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The CP conserved non leptonic K — > it-kit decays are discussed within the chiral quark 
model, including chiral perturbation theory corrections. For the chiral loop correction, 
they arc presented in a way to be more manageable and easier for use. Furthermore a new 
identity for the imaginary part of the chiral loop corrections is derived. All amplitudes arc 
parameterized in terms of quark and gluon condensate and the constituent quark mass. 
The same values for these parameters that can be obtained by a fit of the AI =1/2 rule 
in K — » 2ir give a reasonably good fit in the K — > 37r case. We compare with the work 
of other groups. 
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1. Introduction 

Kaon physics provides us with many interesting phenomena, among them the strik- 
ing AI — 1/2 selection rule in the K — » 777777 decay. The decay of a Kaon into 
three pions has a long history and many attempts have been made to estimate 
its decay amplitudes. The first estimates were done using current algebra meth- 
ods or tree level lagrangian, see Refs. QHHl Afterwards, Chiral Perturbation Theory 
(CHPT) EHH1 a t tree level was used in Refs. EH3 

The first one-loop calculation of the K — > 777777 decay amplitudes in (CHPT) 
was presented in Ref. 1101 but unfortunately the details of this work were lost as 
mentioned in Ref. El Recently, the first full published result, to our knowledge, 
appeared in Ref. 1111 and was later confirmed in Ref. 12 More recently, this was 
followed by including the isospin breaking effects (see Refs. H^HT?>I and references 
therein) . 

The aim of the present work is to give an estimate for K — > 777777 decay amplitudes 
based on the chiral quark mode l (yQM ) 1 1 1 1 approach. The model has been 
applied to the K — > 7777 decays G3H23 an d the K° - K° mixing GEEE1 leading to 
satisfactory results. The strategy is to start with the quark effective Lagrangian, 
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which it turns out splits the physics contributions into short- and long-distances 
encoded in the Wilson coefficients and hadronic matrix elements respectively. 

In evaluating hadronic matrix elements we exploit the chiral quark model as an 
effective link between QCD and chiral perturbation theory, and thus the hadronic 
matrix elements can be estimated using chiral perturbation techniques. The long 
distance contributions computed this way are eventually matched to the short- 
distance Wilson coefficients and the results are compared with the experimental 
values. 

The chiral quark model furnishes us with three input parameters namely the 
constituent quark mass M, the gluon condensate (^f-GG) and the quark condensate 
(qq). A detailed discussion about these parameters and their determination can be 
found in Refs. 123101 

In order for the present paper to be as self contained as possible, we have included 
in the following first two sections the relevant Lagrangian and a brief introduction 
to the xQM. After then the paper is organized as follows: 

in section 4 we present the lowest order chiral lagrangian, while in section 5 the 
isospin decomposition of the decay amplitudes is shown. Section 6 is devoted for 
computing the leading and next-to-leading order for the decay amplitudes. In sec- 
tion 7 we discuss the relevant input parameters for the chiral quark model and in 
section 8 we do the fit with the experimental data and compare with the work of 
other groups. Finally section 9 is for conclusion. We collect all useful formulas and 
relations relevant for our work in the appendix. 



2. The Quark Effective Lagrangian 

The AS = 1 quark effective Lagrangian at a scale /i < m c can be written as 

£ A s=i = — -7=V ud V* s + QiG") • (1) 

v i 

The Qi are effective four-quark operators obtained by integrating out, in the 
standard model, the vector bosons and the heavy quarks t, b and c. A convenient 
and by now standard basis includes the following ten quark operators: 

Qi = (s Q u/3) v _ A (Uf3d a )y_ A , 

Q3.5 = (Sd)y-A T,q (9<?)v T A . / 2 n 
Qi,6 = (s a df3) v _ A Y /q (q0a)vTA , 



My-A T, q e 9 (??) V ±A 



W,9 — 

?8,10 = § (X*gW v _ a J2q e 9 (9/39a)v±A , 



where a, /3 denote color indices (a, /3 = 1, . . . , N c ) and e q are quark charges. Color 
indices for the color singlet operators are omitted. The subscripts (V ± A) refer 
to 7^(1 ± 75) combinations. We recall that Qi^ stand for the VT-induced current- 
current operators, Q3-S for the QCD penguin operators and Q7-10 for the elec- 
troweak penguin (and box) ones. 
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Even though not all the operators in Eq. J2J) are independent, this basis is of 
particular interest for the present numerical analysis becaus e it i s th at employed for 
the calculation of the Wilson coefficients to the NLO in a s . IsSsil-EBI 

In the present paper we will mainly discuss the features related to the first six op- 
erators in Eq. because the electroweak penguins Q7-10 have their contributions 
suppressed by the smallness of their CP conserving Wilson coefficients. 

The functions Zi(fx) and are the Wilson coefficients and Vij the Kobaya- 

shi-Maskawa (KM) matrix elements; r = —VtdV t * s /V u dV* s . The numerical values of 
the Wilson coefficients at a given scale depend on a s , for which we use the recent 
average^a s (mz) = 0.1189±0.002, corresponding to Aq^d = 340±40. We match 
the Wilson coefficients at the mw scale with the full electroweak theory by using 
the LO MS running top mass nit(rnw) — 177 ± 7 GeV which corresponds to the 
pole mass mf°' e = 175 ± 6 GeV. For the remaining quark thresholds we take 
mb(mb) = 4.4 GeV and m c {m c ) = 1.4 GeV. 



3. The Chiral Quark Model 

In order to evaluate the bosonization of the quark operators in Eq. we exploit 
the xQM approach which provides an effective link between QCD and chiral per- 
turbation theory. 

The xQM can be thought of as the mean field approximation of the extended 
Nambu-Jona-Lasinio (ENJL) model for low-energy QCD. A detailed discussion of 
the ENJL model and its relationship with QCD — as well as with the xQM — can 
be found, for instance, in Refs. I41H43I 

In the xQM, the light (constituent) quarks are coupled to the Goldstone mesons 
by the term 



where q T 



-M(q R Zq L + q L S t q R ) , 
(it, d, s) is the quark flavor triplet, and the 3x3 matrix 

2i 



y = exp ( j U(x) 



(3) 



(4) 



contains the pseudoscalar meson octet 



^*0(x) + ^, s (*) 

7T (x) 

K-(x) 



n+(x) K+(x) 



(5) 



The scale / is identified at the tree level with the pion decay constant 

The xQM has been discussed in several works over the years. 111™*! We opted 
for the somewhat more restrictive definition suggested in Ref. 22 (and there referred 
to as the QCD effective action model) in which the meson degrees of freedom do 
not propagate in the original Lagrangian. 



1, 2008 21:52 WSPC/INSTRUCTION FILE threepif 



4 E. I. Lashin 

The QCD gluonic fields are considered as integrated out down to the chiral 
breaking scale A x , here acting as an infrared cut-off. The effect of the remaining low- 
frequency modes are assumed to be well-represented by gluonic vacuum condensates, 
the leading contribution coming from 

<^GG) . (6) 

The constituent quarks are taken to be propagating in the fixed background of the 
soft gluons. This defines an effective QCD Lagrangian £q^ d (A x ), whose propagat- 
ing fields are the it, d, s quarks. The complete xQM Lagrangian is therefore given 

by 

AcQM = £q ff CD (A x ) + ££q M • (7) 

The £ x qm interpolates between the chiral breaking scale A x and M (the con- 
stituent quark mass). The three light quarks (it, d, s) are the only dynamical 
degrees of freedom present within this range. The total quark masses are given by 
M + m q where m q is the current quark mass in the QCD Lagrangian. The Goldstonc 
bosons and the soft QCD gluons are taken in our approach as external fields. A ki- 
netic term for the mesons, as well as the complete chiral Lagrangian, is generated 
and determined by integrating out the constituent quark degrees of freedom of the 
model. The AS — 1 weak chiral Lagrangian thus becomes the effective theory of 
the xQM below the constituent quark mass scale M. In the matching process, the 
many coefficients of the chiral Lagrangian are determined — to the order 0(a s N c ) 
in our computation — in terms of M , the quark and gluon condensates. We neglect 
heavier scalar, vector and axial meson multiplets. 

In conventional chiral perturbation theory the scale dependence of meson loops 
renormalization is canceled out, by construction, through the 0(p 4 ) counterterms 
in the chiral Lagrangian. While in our approach this is maintained for the strong 
sector of the chiral Lagrangian. The tree-level counterterms of the weak sector are 
taken to be fx independent and a scale dependence is introduced in the hadronic 
matrix elements via the meson loops, evaluated in dimensional regularization with 
minimal subtraction. This scale dependence is eventually matched with that of the 
Wilson coefficients. The stability of the final results against fj, variation is numeri- 
cally checked. 

4. Lowest order chiral Lagrangian 

At the leading 0(p 2 ) order in chiral expansion, the strong interaction between the 
SU(3) Goldstone bosons is described by the following effective Lagrangian: 

4?L 9 = Y Tr (D„EtD"E) + y5 Tr (Mtf + EA^) (8) 

where M. is the mass matrix of the three light quarks (u, d and s) and B is defined 
by {<liQj) = —f 2 Bo5ij. The 3x3 matrix E is the same as defined in Eqs. 
The scale / is identified with pion decay constant f n in lowest order. 
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To the same or der, the complete AS = 1 weak chiral Lagrangian in Eq. (JTJ 
can be found in Ref. . Here we just quote the result for the operators Qi-6, which 
turn out to be relevant for the CP conserving amplitude K — > txtxtx. The electroweak 
penguins Qt~\q give a negligible contribution due to the smallness of their Wilson 
coefficients. 

The bosonization of the relevant operators leads to 



(2) 

AS=1 



Gs(Q 3 - 6 ) Tr {XlD^D^Y) + 
G a LL {Qi,2) Tr (A?StD M S)Tr (A^E^E) 
GL(<9i,2) Tr (AiEtD^E) Tr (AjE^D M E) 



(9) 



where A* are combinations of Gell-Mann 5/7(3) matrices defined by (A*-)ft = fe^fe 
and E is defined in Eq. J3J. The covariant derivatives in Eq. © are taken with 
respect to the external gauge fields whenever they are present. 



Table 1. Values of the relevant AS = 1 weak chiral coefficients for two different 
regularization schemes: HV and NDR. The inclusion of the Wilson coefficients 
of the effective quark operators Qi is understood. 



HV 

Gl L {Qi) = -±ti{l~ 



NDR 

G a LL {Qi) = -ikti{i-< 



-f* 

GlliQi) = -ft 
G b 1L (Q2) ■■ 



G a LL (Q2) 

tb ( 
T LL\ 



Hgg)) 



Hgg)) 



G8(Q3) = /^(l-<5 (GG >) 
Gg(Q4) = ft 



Gs(Q 6 ) = 2^fZ (16 



(GG)) 



G a LL (Q2) 

T LL \ 

?L(Q2) = -i-/*(i 



GlliQi) = -f* 
G h TT XQ2) 



°(GG)) 



Gs(Q s ) = f, 



4 1 



'(GG) 



Gs(Q4) = /* 1 - 6 



Gs(Q 5 ) = £ ( 1- - " 



The notation for the chiral coefficients a Gg(<33-6) ; G1 L (Qi^) and G b LL (Qi^) 
reminds us of their chiral properties: G% represents the (8^ x 1 R ) part of the inter- 
action induced in QCD by the gluonic penguins, while the two terms proportional 
to G% L and G b LL are admixture of the (27 L x 1 R ) and the (8 L x 1 R ) part of the in- 
teraction, induced by left-handed current-current operators. These coefficients have 
been evaluated in two different schemes of regularization HV and NDR, and the 
results are given in Table QJ. The parameter <5(gg) is defined as follows, 

N ( a > GG ) 

<W> = y rn^JI' (10) 

where N c is the number of colors. The chiral symmetry breaking scale is identified 
with A x = 2nJfu. 



a It should be understood that the chiral coefficients arc multiplied by the appropriate factors 
containing Fermi constant and Wilson coefficients, but for convenience these factors are not shown. 
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In Table J2J) we report the Wilson coefficients of the first six operators at the 
scales fi = 0.8, \i = 1 GeV in the NDR and HV 75-schemes, respectively. Since 
Rer in Eq. is of O(10 -3 ), the CP conserving part of if — > 7T7r7r transition is 
controlled by the coefficients Zi, which do not depend on m^. 



Table 2. NLO Wilson coefficients at /J, = 0.8 GeV in the NDR and in the HV scheme 
(a = 1/128). The corresponding values at fi = m-w are given in parenthesis. In 
addition one has Z3—g(m c ) = 0. The coefficients Zi(/i) do not depend on mt- 



A (4) 


300 MeV 


340 MeV 


380 MeV 




0.113 


0.119 


0.125 


NDR 


21 


(0.0503) 


-0.591 


(0.0533) 


-0.649 


(0.0557) 


-0.707 


22 


(0.982) 


1.336 


(0.981) 


1.377 


(0.980) 


1.421 


23 


0.0250 


0.0333 


0.0457 


24 




-0.0625 




-0.0799 




-0.1045 


25 


0.0088 


0.0081 


0.0054 


26 




-0.0672 




-0.0881 




-0.1193 


HV 


21 


(0.0320) 


-0.769 


(0.0339) 


-0.879 


(0.0355) 


-1.016 


22 


(0.988) 


1.468 


(0.987) 


1.554 


(0.987) 


1.666 


23 


0.0190 


0.0276 


0.0392 


24 




-0.0392 




-0.0505 




-0.0664 


25 


0.0084 


0.0097 


0.0109 


26 




-0.0372 




-0.0481 




-0.0637 



In comparing with other it is convenient to write the AS = 1 

weak chiral Lagrangian in the following form 

£i 2 Li=28Tr (A^E^E) 



-927 



Tr (A^Z^S) Tr (A^E^E) + - Tr (A^I^E) Tr (A^E^E) 



(11) 



which involves two couplings g$ and <? 2 7 representing respectively the octet (8l, Ir) 
and the twenty-seven (27 l, li?) couplings. There is only one convention used for the 
octet part, however there are different ones for the twenty-seven part, so attention 
should be paid in comparing results. 

In the xQM, these couplings can be written as 

g s _ = ^j=V ud V: s l - {-SGl^QJz^) - 3G2 L (Q 2 )z 2 (,i) + 2G h LL (Q 1 )z 1 ( i x)+ 

2G b LL (Q 2 )z 2 ( f i) + 5 [Gb{Q3)z 3 (jj) + G8(g 4 )2 4 (M) + Gs(Q 5 )z 5 (fj.) + Gs(Q 6 )z 6 (fi)] } 

921 = ^ V "*K,l {(GldQi) + G b LL (Qi)) z 1 ( f x) + (G a LL (Q 2 ) + G b LL (Q 2 )) z 2 ( M )} , (1! 

where G a , Gf, and Gg are the chiral coefficients displayed in Table (|TJ, while Zi(fi) 
are the Wilson coefficients given in Table The coupling gs and g^j_ depend 
implicitly on the parameters (^r-GG), (qq) and M. 
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It seems redundant to express g$ and 527 in terms of three input parameters, 
but one should keep in mind that the complete AS = 1 weak chiral Lagrangian 
contains additional terms transforming as (8^,8^) with their associated couplings. 
Such terms plays an important role in the study of CP violation which is not the 
subject of our present work, a detailed discussion for these terms can be found in 
Ref.EH 



5. K — > tvtvtv Formalism 

There are four distinct channels for K —> tit: it decays: 

K+ -> tt+tt+tt-, (1=1,2), 
K+ -> 7r°7r°7r+, (1=1,2), 

K° -» TT+TT-TT , (7 = 0,1,2), (13) 

K° -> ^°7r 7r°, (I = 1). 



Near each channel we have indicated the final state isospin assuming AI < I. 

In order to write the transition amplitude, it is convenient to introduce the 
following kincmatical variables : 

Si = (pk ~ Pi) 2 and s = -(si + s 2 + s 3 ) = -m|- + to 2 , (14) 

where px and denote Kaon and tt^ momenta (^3 indicates the odd pion i.e the 
third pion in each of the final states 7r + 7r + 7r _ , 7r°7r°7r + , ir + Tr~ir° ). 
We define the dimensionless Dalitz plot variables. 

Y= S -^ and X= S -^. (15) 

Expanding the four amplitudes in Eq. <|13|) in powers of X and Y up to quadratic 
terms, we get 

A++- = A(K + -> tt+tt+tt") = 2a c + (b c + b 2 )Y 

+2c c (Y 2 + X 2 /3) + (4 + d 2 )(Y 2 - X 2 /3), 
A 0Q+ = A{K+ -> TT Tr°Tr+) = a c - (b c - b 2 )Y 

+c c {Y 2 + X 2 /3) - (4 - d 2 ){Y 2 - X 2 /3), 

A+- Q = V2A{K° -» tt+tt-tt ) = a n - 6„F - h 2 X 

+c n (Y 2 + X 2 /3) - d n (Y 2 - X 2 /3) + U 2 XY, 
Aqoo = V2A(K° - TT Tr 7r°) = 3 [a„ + c„(F 2 + X 2 /3)] . (16) 



Here, we follow the notations and conventions of Refs. 

A representation al terna tive to that in Eq. 1)160. which has been adopted in fits 
to experimental data, E33I i s the expansion in terms of amplitudes with definite 
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isospin selection rules. Such an expansion can be written as 



A(K L 




(ai + a 3 ) - + fo)Y 

+ (Ci - 2C 3 )(F 2 + X 2 /3) + (d - 2£ 3 )(F 2 - X 2 /3), 



A(K L 
A(K+ - 




3(ai + a 3 ) + 3(d - 2( 3 )(Y 2 + X 2 /3) 7 
-(2a, - a 3 ) - (A - i& + %/3 73 )F 



7T + 7T + 7r ) 



+ 





-2(Ci + ( 3 )(Y 2 + X 2 /3) + (6 + & " C 3 )(^ 2 
-i(2a a - a 3 ) + (#l - ^ 3 - V3 73 )F 

-(a + c 3 )(r 2 + x 2 /3) - (a + e 3 + e 3 )(^ 2 - 

2 a- 4 , 

- V3 l3 X - -&XY. 



X 2 /3), 



X 2 /3) 



(17) 



In Eq. I|17p. the subscripts 1 and 3 refer to AI = 1/2 and AI = 3/2, respectively. 
The relation between the amplitudes in Eq. I|16|) an d those in Eq. (|17fl is easily 
found to be 



6. K — > tvtvtv Amplitudes. 
6.1. Leading order 

Employing the lowest order chiral Lagrangian given by Eq. © and Eq. ©, we can 
calculate the K — > nmr amplitudes. The number of the contributing diagrams turn 
out to be huge, more than three hundred diagrams. The complete relevant Feynman 
rules are too lengthy to be mentioned but it can be found in Ref. 03 Necessary 
integral formulas for calculating Feynman graphs arc collected in appendix A along 
with many other useful formulas. 

The contributing diagrams fall into two classes : 

• (i) Tree diagrams as given in Fig. Q turn out to be easily calculated and 
their results are polynomial functions of the external momenta. 

• (ii) Loop diagrams, The ones in Fig. (J2J are tadpole-kind corrections and 
their results are polynomial functions of the external momenta. The rest of 
the diagrams in Fig. with at most two insertion of strong chiral Hamil- 
tonian are the most complicated ones. Their results are non polynomial 
functions of X and Y variables defined in Eq. H15J1 . Hence, we expand the 
results in Taylor series around the point (X = 0, Y = 0) which is the center 
of the Diltaz plot variables. 



d c 
hi 



-a,\ + a 3 /2, a n 

ft/2, b n 

—& - Ca, c n 

Ci+6) d n 

-V3~7 3 , d 2 




(18) 
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Fig. 2. One-loop diagrams of tadpole kind correction for K — > tttttt. The black box and circle 
indicate the insertion of the weak AS = 1 and strong chiral Hamiltonian respectively. 
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Fig. 3. One-loop diagrams with horizontal loop for K — > nnn. The black box and circle indicate 
the insertion of the weak AS = 1 and strong chiral Hamiltonian respectively. 



The tree diagrams in Fig. (1) give the following results. 



A ++ -(Ga) 

Aoo+ (G a ) 

A + -o(Ga) 

Aaoo(G a ) 



— A (m 2 K + 3mlY) , 

1 (2m 4 K - <om\Y + m 2 K m\ (-2 + 15Y)) 
6 f 4 m 2 K — m 2 

1 2 (-Sm 2 K (X + Y) + 2ml {X + 3Y)) 



2/4 
0, 



rn 



m" K -m* 



(19) 



A ++ -{G b ) = ^mlY, 



Aoo+(Gb) 
A+-o(G b ) 



3 m 2 K mlY 
2 f 4 m 2 K — m 2 ' 

1 2m 4 K + 6m 4 (X + 2Y) - m 2 K ml (2 + 9X + 3Y) 
6 f 4 m 2 K ~ml 



A 00 o(G b ) = , (20) 
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1 

37 4 

2 



A ++ _(G 8 ) = — (-2m! + 3m^y 
1 



A 00+ (G 8 ) = -jjz (m 2 K + imlY) , 
A + . (Gs) = ^(m% + 3mlY), 

Aooo(G 8 ) = . (21) 

Expressing our leading results in terms of amplitudes with definite isospin se- 
lection rule defined in Eq. (|17|l . we find for the AI = | amplitudes 

ai = J 1 ^ { 9 ~ ~ \ 9 ^) ' Pl = J 1 ^~ m ^ ( 5 - " \ 52Z ) ' (22) 



and the following for the AI - " 



2 

1 10 o „ I m'i ,_ ., „ 5 

2~ 

K 



f*Z7"' KySb ^~T* ml -m* VJ "' K ""^18 



327, /3 3 = —32 — ZL 3T ( 5to K ^ 14rr 4) ^ ff2 7, 

^ = ^^a(3-|-2^)-=|^ (23) 

where the masses m ff , rn^ and m, are taken to be 0.138, 0.498 and 0.548 GeV 
respectively. 

The leading order results in Eqs. H22J1 -123 JI agree with those found in Refs. 1101 
1111 1441 after taking care of the used different conventions. 



6.2. Next-to-leading order 

In our present work, only one-loop correction is included, the other next-to-leading 
contribution coming from the order 0(p ) counter terms of weak chiral Lagrangian 
will not be considered here. In fact, these contributions of weak counter terms can 
be computed using the xQM which is a rather lengthy calculation. We leave this to 
a future work. 

Other kinds of one-loop corrections are not shown graphically in Figs. i|2!3|) 
namely, wave function and f renormalization. Such kind of corrections originate 
purely from the strong sector. 

The wave-function renormalization which arise in the chiral perturbation from 
the direct calculation of K — > K and tt — > tt propagators are given at 0(p 2 ) by 



Z k = 1 + i 
Z n = 1 + i 



-h{m^) + -J 2 K) 



2 1 1 

-h{m n ) + -h{m k ) + -h{m K ) 



(24) 
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where ^im) is 



h(m) 



16tt 2 



[l-log(m 2 )+log(/! 2 )] 



(25) 



To include the wave-function renormalization for the amplitude K — > tttttt, each 

i i 

tree level result should be multiplied by the factor Z£ Zj^. 

Concerning the / renormalization; at the lowest order in chiral perturbation 
there is no distinction between the decay constants. After considering one-loop 
correction, the physical decay constants are identified according to the following 
equation 



f K = / 1 



3 3 3 

-h{m^) + -h{m K ) + g-fe(m^) 



P 

U = J < i 



him*) + ^h{m K ) 



(26) 



where I2 (m) is defined in Eq. I|25l) . The inclusion of the / renormalization is made by 
expressing l// 4 appearing at the tree level in terms of the physical decay constants 
U(= 0.0924 GeV) and f K (= 0.1130 GeV). One power of / is eliminated in favor of 
/k while the others in favor of f v through using the relations in Eq. (|26J) . 

In regularizing our loop integrals, we used the dimensional regularization with a 
modified minimal subtraction scheme (MS) h . For the sake of convenience, a factor 
10~ 6 is extracted from the results. A part from the tree level result, the one-loop 



b In this subtraction scheme all the combination 
Euler constant. 



7b + log 4-7T is subtracted, where is the 
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results after including wave function and / renormalization are 



A ++ _(G a ) = [(2835. + 183.3 1+ 1112. log^ 2 ) + 

Y (471.1 - 107.4^ + 217. log^ 2 ) + 
X 2 (4.259+ 1.317^+ 1.053 log M 2 ) + 
Y 2 (-15.02 - 6.924i - 5.455 log M 2 )] , 

A 00 +(G a ) = [(1418. + 91.66 1 + +555.9 log M 2 ) + 

Y (733.6 + 71.1 i + 290.6 log^ 2 ) + 
X 2 (4.352 + 2.937 i + 0.7177 log^ 2 ) + 
Y 2 (-14.18- 10.3 i- 3.301 log M 2 )] , 

A+-o(G a ) = [(648.6 + 245.5 log/i 2 ) + 

X (-401.5 + 12.1 i - 169.2 log fj 2 ) + 

Y (-423.7- 10.79 i- 163.6 log^ 2 ) + 
X 2 (-6.657- 1.486 i- 2.632 log^ 2 ) + 
Y 2 (-3.824 - 6.541 i - 0.7177 log^ 2 ) + 
X Y (-18.35 - 9.995 i - 4.689 log fi 2 )] , 

Aooo(G a ) = [(1946. + 736.5 log M 2 ) + 

X 2 (-11.9- 5.499 i- 4.306 log^ 2 ) + 
Y 2 (-35.69 - 16.5 i - 12.92 log M 2 )] , 



(27) 
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A ++ _(G b ) = -L [(1028. + 441.1 log/i 2 ) + 

Y (771.4- 12.75 i + 310.5 log/i 2 ) + 
X 2 (1.764 + 0.492 i + 0.1914 log fj 2 ) + 
Y 2 (-29.85 - 12.47 i - 9.761 log fj 2 )] , 

A 0Q+ (G b ) = [(513.8 + 220.6 log M 2 ) + 

Y (432 .4 - 23.54 i + 196.7 log /i 2 ) + 
X 2 (1.27 + 1.756 i - 0.1435 log fi 2 ) + 
Y 2 (-16.09- 10.77 i- 4.163 log fj 2 )] , 

^+-o(G b ) = [(1555. + 91.66* + 581.9 log M 2 ) + 

X (-401.3 +12.H- 169.1 log M 2 ) + 

Y (-123.1 + 83.85^- 69.95 log/i 2 ) + 
X 2 (-3.575 - 0.3046 i- 1.77 log fi 2 ) + 
F 2 (-1.916- 6.072 i + 0.1435 log^ 2 ) + 
XY (-18.35- 9.995 i- 4.689 log^ 2 )] , 

A oao (G h ) = [(4665. + 275. i + 1746. log M 2 ) + 

X 2 (-6.321 - 3.493 i- 2.584 log/i 2 ) + 
Y 2 (-18.96- 10.48 i- 7.751 lo gy u 2 )] , 
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A ++ _(G 8 ) = i [(-1808. - 183.3 i- 670.6 log/i 2 ) + 

Y (300.3 + 94.64 % + 93.49 log /t 2 ) + 

X 2 (-2.495 - 0.8252i - 0.8612 log// 2 ) + 
Y 2 (-14.82 - 5.551 % - 4.306 log// 2 )] , 

,4 00+ (G 8 ) = i [(-903.9 - 91.66 i- 335.3 log/?) + 

y (-301.2 - 94.64 i- 93.85 log// 2 ) + 
X 2 (-3.082- 1.181 i- 0.8612 log// 2 ) + 
F 2 (-1.908- 0.4689 i- 0.8612 log/i 2 )] , 

A+_ (G 8 ) = i [(906.6 + 91.66 i + 336.4 log/?) + 

X (0.2895 + 0.1209 log/i 2 ) + 

Y (300.6 + 94.64 i + 93.61 log /t 2 ) + 
X 2 (3.082 + 1.181 i + 0.8612 log/i 2 ) + 
Y 2 (1.908 + 0.4689 i + 0.8612 log /t 2 )] , 

A)oo(G 8 ) - [(2720. + 275. i + 1009. log// 2 ) + 

X 2 (5.577+ 2.007 i+ 1.722 log// 2 ) + 

y 2 (16.73 + 6.02 i + 5.167 log// 2 )] , (29) 

where fx is the renomalization scale measured in units ol GeV and i = \f— 1. 

As a check of the calculations, one can find different ways for determining the 
parameters a c , b c , ■ ■ ■ etc. defined in Eq. Ijl6|l which should give the same result. As 

an example a c can be determined from both and Aqq + leading to the same 

result. Similar patterns hold for the rest of the parameters. 

Another check comes from observing that the term 

Tr (A^E^E) Tr (A^Z^E) - Tr (A^Z^E) Tr (A^Z^E) , (30) 

transform as an octet (8i x 1#). Computing Feynman diagrams for each part sepa- 
rately and then subtracting them, the result should be the same as of the octet part 
of the weak chiral Lagrangian. In fact, these two tests constitute non trivial checks 
for the performed calculations and our calculations did pass them successfully. There 
is an additional check which holds for the imaginary parts of the amplitudes and is 
explained in appendix A. 

Expressing our one-loop results in terms of dimensionless amplitudes with def- 
inite isospin selection rules (apart from the tree level results which are given by 
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Eqs. I|25 )l -ip5 )l . we get for the octet 



(8) 


= 12.402- 


f 1.256 i + 4.601 \og(fi 2 ) 




= 0.025 + 


0.010 log(/i 2 ) , 


/#> 


= -4.122 


- 1.297 i- 1.284 log(> 2 ) 


M 8) 


= 0.001 + 


0.00055 log(> 2 ) , 


C (8) 


= 0.076 + 


0.028 i + 0.024 log(^ 2 ) , 




= -0.050 


-i 0.021 -0.012 log {n 2 ) 


7s W 


= 0.003 + 


0.0014 log(> 2 ) , 



and for the 27-plet 



(27) 


= -4.250 


- 0.134 i 


- 1.995 log(/x 2 ) , 


4 27) 


= 31.494- 


\- 1.396 i - 


f 12.216 log(/Lt 2 ) 




= 1.104 + 


0.144i + 


0.522 log(/i 2 ) , 


/f 7) 


= 4.455 - 


1.195? + 


1.932 log(/i 2 ) , 




= -7.941 


+ 0.239 i 


- 3.346 log(/i 2 ) , 


C (27) 


= -0.006 


-i 0.003 


- 0.0004 log (/i 2 ) 


a(27) 
S3 


= 0.095 + 


i 0.048 + 


0.037 log( M 2 ) , 


^(27) 


= 0.009 + 


i 0.002 + 


0.003 log(^ 2 ) , 


t(27) 
S3 


= -0.063 


+ i 0.024 


- 0.033 log (/i 2 ) 


t'(27) 
S3 


= 0.314 + 


i 0.171 + 


0.080 log (> 2 ). 



The results in Eqs. (|31[) - l|32(l are made to be dimensionless through multiplying by 
the factor / 2 . Then the full one-loop contributions to the isospin amplitudes are 
obtained by setting = (f^j-J + (ff)^- 

The values of the isospin amplitudes depend on the scheme and scale of renor- 
malization. For a meaningful comparison we restrict ourselves to the imaginary part 
which are scale and subtraction-scheme independent, leaving the comparison of the 
real parts for the final fit to the experimental results. 

The results in Table indicate that there is a reasonable agreement of our 
calculations and those of Ref. for most of the imaginary parts of the isospin 
amplitudes. The rest of the parameters p[ 8 \ and 73 ' disagree. 

In Ref. 1481 the imaginary part of the amplitude of the channel K + — > Tr + ir + ir~ 
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Table 3. Comparison of the imaginary parts of the isospin amplitudes. 





(8) 


M 8) 


A8) 
SI 


5 (8) 




/3f 7) 


SI 


? (27) 


Kambor et al. 


1.00 


0.50 


0.019 


0.00 


- 0.67 


0.33 


- 0.012 





Our calculations 


1.26 


-1.30 


0.028 


- 0.021 


-0.42 


0.432 


-0.009 


0.006 






/f 7) 


(27) 
% 


a(27) 
S3 


.(27) 
S3 


,,(27) 
S3 






Kambor et al. 


6.70 


-2.19 


1.12 


0.066 





0.43 






Our calculation 


4.2 


- 3.585 


0.717 


0.144 


0.072 


0.513 







was calculated for the octet part of the weak chiral Lagrangian leading to the result 

Im [A (K+ -> TT+TT+TT-)] = 



-1 



192 7T /« 



/ 4?n^ 

6 4/1 — (m 2 K + ml -si- 8 2 ) (2ml - *i - s 2 ) 



S3 



/ Am?* 

+2 Wl * (9m£ + (2rn 2 + Sl ) + Sl ( 5si + S2 ) _ m 2 ( llsi + 4fi2 )) 



8 1 



/ 4 771 ^ 

2 W 1 - — £ (9mJ + TO |_ ^ 2m 2 + s ^ + S2 (5s2 + Si) _ m 2 (1Lsa + 4si) ^ 



(33) 



Expanding the result in Eq. (|33|l in terms of X and Y up to quadratic terms, we 
get 

Im [flA{K+ -> tt+tt+O] = -2. 51 + 1. 2973^-0.0113109 X 2 - 0.0760878 y 2 . (34) 

The imaginary parts of the dimcnsionlcss isospin amplitudes determined from 
Eq. (jSU are 

af ] = 1.26, (3f } = -1.30, c[ 8) - 0.028, f< 8) = -0.021, (35) 

which are in a full agreement with our result in Table J3J . A similar agreement holds 
for the twenty-seven part after paying some caution to the convention followed 
in Ref. 021 Our analytical result of the Im [A (K + — ► ir + ir + Tr~)] , for both octet 
and twenty seven parts of the Lagrangian, with those of Ref. 1481 can be found in 
appendix A. 

It is difficult to point out the source of discrepancies between the two results 
mentioned above, in particular, since the details of work Ref. E3 were lost as stated 
before in the introduction. However, as shown in the appendix, some theoretical 
identities relating the imaginary parts of the amplitudes are satisfied for our results 
but fail in Ref. 1101 Moreover, our results agree with other previously calculated 
ones in Ref. 0HI Thus, we suspect that the discrepancies might be attributed to 
some missing diagrams in the calculations of Ref. E3 
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7. Input parameters 

Before discussing the fit of the isospin amplitudes with the experimental data, we 
should consider the relevant input parameters for the chiral quark model. 

The value of the input parameters M, the gluon condensate (^GG) and the 
quark condensate (qq) are chosen such that to fit the AI = 1/2 selection rule 
characterizing the I — and 1 — 2 amplitude of the non-leptonic Kaon decay 
(K -► 2 vr). 

Such a fit was done in Ref.l2f>lfor the complete 0(p 4 ) calculation in the framework 
of xQM. Since, in our present work, we stop at 0(p 2 ) order while including the 
meson loop correction, we redid the fit in order to determine the appropriate values 
of the parameters. The resulting values of the input parameters which fit the AI = 
1/2 rule up to 20% are 

( — GG) = (358^ 2 MeV) 4 , 

(qq) = -(260±« o MeV) 3 , 
M = 200 ± 20 MeV. (36) 

Let us focus on the coupling constants gs and c/27 defined in Eq. 1)12(1 . For 
convenience we divide them by /~ to produce dimensionless quantities measured 
in units of 10~ 8 . For our central values of the parameters (M, (^-GG), (qq)), the 
couplings take the following values 

|§ = 3.42 , = -0.382. (37) 

J TV J 7T 

These values show clearly the dominance of the octet coupling (gs), which is neces- 
sary for producing the A I = | selection rule. 

We illustrate the dependence of the couplings on the input parameters in Fig. J2J. 
For convenience we kept (qq) and M fixed at their central values in Eq. I|36|l . while 
we varied the gluon condensate in the following range 

(-^GG) = (329 - 423 MeV) 4 , (38) 

As evident from Fig. increasing the gluon condensate tends to enhance the 
coupling g$ while suppressing the coupling (727- It is only the coupling gs which 
depends on the quark condensate and M. It turns out that the coupling <? 8 is 
enhanced by increasing the quark condensate or by reducing M. In the foregoing 
discussion, the renormalization scale (i was fixed at 0.8 GeV, and the scheme is 
taken to be the HV one. The 75-scheme and scale dependence of our final results 
are discussed at the end of the next section. 

It is worth mentioning that the phenomenological determination of the quark 
and gluon condensate are a complicated issue and the literature offers different 
determinations of them. To fix the ideas, the condensate entering our computation 
are identified with those obtained by fitting the experimental data using the QCD 
sum rules (QCD-SR) methods or the lattice computations. 
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5^ — , — ■ — ■ — ■ — , — ■ — ■ — ■ — , — ■ — ■ — ■ — , — ■ — ■ — ■ — ^ 



-g27*10 



0.34 



0.36 0.35 

<GG> 



0.4 



. 42 



Fig. 4. The solid (dotted) line represents the coupling gg/f% (1 927//^ I X 10) as a function of the 
gluon condensate, (qq) and M are fixed at their central values in Eq. 1361 . 

Regarding the gluon condensate, it was first determined by Shifman, Vainshtein, 
and Zakkharov fSVZ) using QCD-SR. Later, lattice QCD was used for this 

determination. LiiHMI We prefer to take the value of the gluon condensate in the 
range 



-GG) = (376 ± 10 MeV) 4 , 



(39) 



which encompasses the result of QCD-SR analysis . ISSHSII The corresponding ranges 
provided by the lattice calculations in Refs. EU EH are respectively, 



[—GG) = (460 ± 21 MeV) 4 , and ( — GG) = (651 ± 40 MeV) 4 

7T 7T 



(40) 



and they are larger than the estimates provided by QCD-SR. 
For the quark condensate, we take the range 

-(200 MeV) 3 < (qq) < -(280 MeV) 3 (41) 

in order to include the central values and errors of the QCD-SRpHl and the lattice 
estimate™!] 

As to the constituent quark mass M, we consider the range 

M w 200 - 250MeV, (42) 

which is consistent with the estimate based on processes involving mesons. E3I Such 
a value is smaller than the value M ss 330 MeV often quoted from baryon physics. 



8. Comparison with the experimental data 

In comparing the K — > 7r7T7r amplitudes with the experimental data, we use the 
expansion given in Eq. Ijl7(l . The analysis of the parameters at, a.2, ■ ■ ■ etc. was first 
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made by Devlin and Dickey in the late seventies.^ In the early nineties the analysis 
has been redone by Kambor et al. .E3 using recent data and including the AI = 3/2 
quadratic slope parameters £3 and £3 in Eq. l(T7)l. More recently ^2 a f u ll fit has been 
done using recent data in Ref. 1641 The fitted values for the parameters cti, ct2, ■ ■ ■ etc. 
are displayed in Table (@J. 



Table 4. Determination of K — > 3n isospin amplitudes, in units 
of 10" 8 . 



Quantity 


Rcf. 3 


RefHOl 


Ref.Sl 


ai 


91.4 ±0.24 


91.71 ±0.32 


93.16 ±0.36 


a s 


-7.14 ±0.36 


-7.36 ±0.47 


-6.72 ± 0.46 


Pi 


-25.83 ± 0.41 


-25.68 ±0.27 


-27.06 ± 0.43 


P8 


-2.48 ± 0.48 


-2.43 ±0.41 


-2.22 ± 0.47 


73 


2.51 ±0.36 


2.26 ± 0.23 


2.95 ±0.32 


Ci 


-0.37 ±0.11 


-0.47 ±0.15 


-0.40 ±0.19 


C3 




-0.21 ±0.08 


-0.09 ± 0.10 


£l 


-1.25 ± 0.12 


-1.51 ±0.30 


-1.83 ± 0.30 


?3 




-0.12 ±0.17 


-0.17 ±0.16 


?3 




-0.21 ±0.51 


-0.56 ± 0.42 



As discussed in Refs. 23: 42~Tl it was shown that within the frame work of xQM , the 
AI = 1/2 rule in the case of K — > 2ir decay can be reproduced. Here we show how 
this can be achieved for the case K — > 37r. Since our results depend on the following 
parameters Gluon (^-GG) , quark (qq) condensate and the constituent mass M, 
we choose the values which fit the AI = 1/2 rule for the K — > 2n determined by 
Eq. 1)36(1 as mentioned in section 7. We show the numerical results for the isospin 
amplitudes at the central value of this range in Table 101 ■ The results are clearly 



Table 5. The isospin amplitudes, in units of 10 at the central value of the range 
given by Eq. 



ai 


03 


01 




73 


Ci 


Cs 


€1 


6 




64.15 


-13.23 


-18.8 


-2.18 


4.02 


0.23 


-0.03 


-0.16 


0.018 


-0.11 



off by a factor of two for 013 while being in a reasonable agreement for the rest of 
the linear slope parameters. 

In order to improve the result, there should be a suppression of 013 and 73 which 
can be achieved by lowering the coupling 327- This can be done by increasing the 
gluon condensate. We find the best fit corresponds to the following values; 

(-GG) = (376 MeV) 4 , 

(qq) = -(260 MeV) 3 , 
M = 200 MeV. (43) 

which for the gluon is slightly out of the range specified by Eq. Q36Jl ■ but it is still 
safe for fit of the A J = 1/2 in the K — > tttt decay within 30%. 
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Here in Table 10 we list our results together with the ones obtained in Ref. ITTI 
corresponding to the case of neglecting higher order counter terms for the sake of 
a meaningful comparison. As it can be seen there is a reasonable agreement within 



Table 6. The value of K -* 3?r 
isospin amplitudes, in units of 10 — 8 
for case of neglecting higher order 
counter term. 



Quantity 


Ref. 11 


Our results 


ai 


59.4 


68.55 


a 3 


-6.5 


-9.91 


Pi 


-21.9 


-20.11 


fti 


-1.0 


-1.64 


73 


2.5 


3.02 


Ci 


0.26 


0.24 


Cs 


-0.01 


-0.02 


6 


-0.46 


-0.17 


£s 


-0.01 


0.01 


?3 


-0.06 


-0.08 



30% for the parameters ct\, (3i, 73, Ci and £3. Even if we compare with the full fit^^ 
listed in Table (J2J, there is a reasonable agreement for the linear slope parameters 
an, 0i, 03, and 73 except for 03 which agree within 50%. As to quadratic slope 
parameters only £3 is compatible with the fit while the others still need further 
improvement. The result are expected to be improved after including the higher 
order counter terms. These kind of correction can be calculated in the chiral quark 
model, which are rather lengthy and will be the subject of future work. Another 
feature which is worthy to be discussed namely the 75-scheme and scale dependence 
of the amplitudes. As a direct measure of the scale dependence we define 



A P A = 2 



A(0.8 GeV) - A(1.0 GeV) 



A(0.8 GeV) + A(1.0 GeV) 



(44) 



while the difference between the HV and NDR results is quantified by 



A 75 Ai 



^NDR _ ^HV 



^NDR 1 ^HV 



(45) 



As evident from Table 10 the scale dependence of the isopin amplitudes is within 
20% excluding £3 and £3 for which the scale dependence turns out to be about 30%. 
As to the 75-scheme dependence is not large for the amplitudes and remains below 
10%. The overall dependence is not large and still under control. 

Finally one can say, with caution, that chiral quark model provides a coherent 
right picture for the K-meson physics not only for K — > irn decay as shown in 
fiefs. !23H27l but also for case K — > 777777. 
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Table 7. The scale and 75 scheme dependence of the isospin amplitudes A (in units of 10 8 ). The amplitudes 
are computed for the values of the parameters given by Eq. 



A 

ai 


At = 0.8 GeV 


H = 0.9 GeV 


H = 1.0 GeV 


A M A 


HV 


A J5 A 


NDR 


HV 


A 75 A 


NDR 


HV 




NDR 


HV 


NDR 


68.55 


4% 


66.05 


61.72 


3% 


59.66 


57.77 


3% 


56.16 


17% 


16% 


«3 


-9.91 


7% 


-10.70 


-11.18 


6% 


-11.91 


-12.30 


6% 


-13.0 


22% 


19% 


ft 


-20.11 


4% 


-19.37 


-18.05 


3% 


-17.44 


-16.84 


3% 


-16.37 


18% 


17% 


ft 


-1.64 


8% 


-1.77 


-1.85 


6% 


-1.97 


-2.03 


5% 


-2.14 


21% 


19% 


73 


3.02 


8% 


3.26 


3.38 


6% 


3.60 


3.69 


6% 


3.9 


20% 


18% 


Ci 


0.245 


4% 


0.235 


0.224 


3% 


0.217 


0.213 


3% 


0.207 


14% 


13% 


Cs 


-0.022 


9% 


-0.024 


-0.026 


7% 


-0.028 


-0.030 


3% 


-0.031 


31% 


25% 


£1 


-0.168 


4% 


-0.162 


-0.151 


3% 


-0.146 


-0.142 


3% 


-0.138 


17% 


16% 


€3 


0.014 


7% 


0.015 


0.017 


6% 


0.018 


0.019 


5% 


0.020 


30% 


29% 


«3 


-0.080 


7% 


-0.086 


-0.089 


5% 


-0.094 


-0.096 


6% 


-0.102 


18% 


17% 



9. Conclusion 

In this paper we have calculated the K — > tttttt amplitudes by including the one 
loop correction in chiral perturbation theory in the framework of chiral quark model. 
Our calculations are independent of what have been carried out in Refs.EB^2 The 
results for the chiral loop corrections are presented in a way to be more manageable 
and easier for use. Further more we derive a useful new identity for the imaginary 
part of the K — > nmr amplitudes, the details are left for the appendix. 

Our results for the isospin amplitudes are compared with the recent fit done 
by, we find a reasonable agreement for the linear slope parameters except a^. 
As to the quadratic slope parameters only (3 is compatible with the recent fit. 
However restricting the comparison to the fit while neglecting higher order counter 
terms leads to a reasonable agreement for almost all the slope parameters. 

The result can be further improved by including higher order counter terms in 
the strong and weak chiral lagrangian, which can be calculated in the frame work 
of chiral quark model. However the calculations are rather lengthy and will be the 
subject of future work. 
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Appendix A. Integral Formulas, Conventions and Imaginary Part 
of the Amplitudes 

In this appendix A we collect all the integral formulas necessary for our calculations 
and summarize the conventions used in writing the weak chiral Lagrangian. More- 
over we present some of our calculated analytical results for the imaginary part of 
the amplitudes and prove some identities which should be satisfied by the results. 
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All the loop integrals are evaluated by using dimensional regularization assuming 
the space- time dimension d = 4 — 2 e. In subtracting the divergences we follow the 
modified minimal subtraction scheme (MS), which means that not only the pole 
term | is subtracted but the whole combination (7 — 7e + log47r). 

The necessary formulas for loop integrals are: 

fi (x, mi , m 2 ) = m\ - x (m\ - to?) , 
f 2 (x,p 2 ,m 1} m 2 ) = m\-x(\ - x)p 2 - x [m\ - to?) , (A.l) 
h (x,y,p 2 ,m 1 ,m 2 ,m 3 ) = m| - x(l - x)p 2 - x [m\ - to?) -y{m\- m 2 ) . 

The finite parts for the required loop integrals are: 
I- ddq 1 " 



^m^l-lnm^+ln^), 

(2^ [( 9 +p) 2 -m2j(,2-m|) = 

jff^z fa dx\nf 2 (i,p 2 ,mi ) m 2 ) +ln/i 2 J 



|^L_ 1 



f d"q (q-Pl) 

(2^F [(q+p) 2 -m 2 ](q 2 -m 2 ) 



J d g (g-Pl)(g-P2) 



[jg 1 cfa x p ■ pi ln/2 (i,p 2 ,mi.™2) - P ■ Pi ln/i 2 J , 

(2^ [(q+p) 2 '-'m\)(q 2 '-ml) = 

32^ Pi 'P2 / dxf 2 (x,p 2 ,m 1 ,m 2 ) (l - ln/ 2 (x,p 2 ,m 1 ,m 2 ) - In ^ 2 ) + 

T ^p-pip-p 2 [j,, 1 ^!: 2 (-ln/ 2 (x,p 2 ,mi,m 2 ) +ln/i 2 )j , (A.2) 
dfg 1 





K<?+p) 2 


-m 2 ](q 2 -m 2 )( 9 2 


-m 2 ) 
16^ J( 


d d 9 




(<3'Pl) 




(2,r) d 


[(9+p) 2 


-m 2 ](q 2 -m 2 )( 9 2 

167T 2 


-m%) 
p-pi 


d d q 




(<?'Pl)(<?'P2) 




(27T) d 


[(■j+p) 2 


-m 2 ](q 2 -m 2 )(g 2 


-m 2 ) 



"1 J fl-I 



/ 

T ^j-p-pip-p 2 /gdi Jq~ x dy x 2 f.3 1 (x,y,p 2 ,m 1 ,m 2 ,m 3 ) 
+ 3^rPl-P2 [-Jo 1(l!:E fo~ X dy ln /3 (i,t/,p 2 ,mi,m2,m 3 ) + ln^ 2 ] . 

It should be understood that the above result for loop integrals are used after 

expanding them properly about the origin of the Diltaz plot variables (X = Y = 0) 

up to quadratic terms in X, Y. 

The integral formulas needed to compute analytically the imaginary part of the 
Feynman diagrams are: 

( 2 ^) 2 / j0*Hl 2 -™l) S[(a + p) 2 -ml] = 



iL v/( 1 ~^)' 



(27ri) 2 / *A q ■ pi S (q 2 - ml ) 5 [(, + p) 2 - m 2 ] 



P-Pi, (A.3) 
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(2niff g . pi q . p2 6(q*-ml) S [(, + p) 2 - ml] 



im l V 2 1 f 1 ml \ I ft 4m2 ' 

192^1--^) V 2 - PI -P2 - W ( 1 - ^ ) ~ -pF) P-P1P-P2, 



where S denotes the Dirac delta function. 

The weak chiral Lagrangian according to the convention followed in Ref@Hl is 

4 2 Li = 9s Tr (XlD^D^E) 

+327 [Tr (A^Et^S) Tr (AfE^E) + Tr (A^E^E) Tr (A^Z^E) 
- Tr (A^E^E) Tr (AlE^E)] . 

(A.5) 

The SU(3) projection for the twenty-seven part are 1^1 

|27, -) = Tr (A^E^S) Tr (AfE^E) + 4 Tr (A^E^E) Tr (A^L^E) 



-5Tr (Af E^E) Tr (Xp^D^E) 
f (A^E^E) Tr (AfE+D^E) + 
Tr (AlEt^E) Tr (A^E^E) 



|27, ~) = Tr (AJjE^E) Tr (AfE+D^E) + Tr (A^E^E) Tr (A^Z^E) 



(A.6) 

Therefore, we have 

l22) = §|2r,§> + i|2 Z ,i>, 

= H Tr (A^Et^E) Tr (A^E^E) + Tr (AfE^E) Tr (AjE^E) (A.7) 

So, only the |27, §) component was considered in Ref . 1481 This component induces 
the AI = 3/2 transition. 

The full imaginary part of the amplitude K + — > 7r + 7r + 7r~ for the octet part of 
the weak chiral Lagrangian given by Eq. 

is 



Im [AW (K+ ^Tr + TT+TT-)] 



192 n /| 



6 W 1 — 4 (m 2 K +m1 - si — s 2 ) (2 ml - s 1 - s 2 ) 

S3 



+2 \j 1 - (9 + (2 m 2 + Sl ) + si (5 si + s 2 ) - m 2 (11 si + 4 s 2 )) 



(A.8) 
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and for the twenty-seven part 

Im[A< 27 > (K+ -^n+n+n-)] = 



192 7T/6 V S3 



" ' 12 Jl - (rn\ + ml - Si - s 2 ) (3 m\ + 8 m 2 - 4 ( Sl + s 2 )) 



' [m| (23 m 2 _ n Sl) 



si 
,2 



+m 2 (31 + 2 m 2 (3 si - 32 s 2 ) - si (si - 16 s 2 )) 
-2 m 2 n (45 + si (13si + 8s 2 )-2mJ (17 si + 16 s 2 ))] 

: "' S [m| (23 m 2 - 11 S2 ) 



"K ~ '"t V S2 



+m£ (31 m* + 2 m 2 (3 s 2 - 32 si) - s 2 (s 2 - 16 si)) 
-2 m 2 (45 m% + s 2 (13s 2 + 8 si) - 2 m 2 (17 s 2 + 16 si))] } . 

(A.9) 

The authors of Ref!^ did neglect terms of O(rnJ) in the weak vertex. But their 
result for the octet parts agrees with ours. For the twenty-seven part, there is 
disagreements in terms proportional m l, w hich can be safely neglected. Here is the 
result found for the twenty seven part. 

Im[A( 27 > (K+ -^n+n+n-)] = 



192 7T/6 V s 3 



" 1 12 J 1 - (mjc + ml - si - s 2 ) (3 m\ + 8 m 2 - 4 (si + s 2 )) 



1 - (-63 ml + s 2 (-16 si + s 2 ) + 2 m 2 (32 si + 7 s 2 ) + ml (-23 m 2 + 11 s 2 )) 

si 



1 - (-63 + si (-16 s 2 + s x ) + 2 m 2 (32 s 2 + 7 si) + m| (-23 m 2 + 11 si)) 



(A.10) 



Some identities for the imaginary part could be derived, but first let us simplify 
our notations as follows 

Tr (A^E^S) Tr (A?Et£>"E) -» G a , 
Tr (AfEtD^E) Tr (AjS+DT) -► G b , 
Tr (A^E^E) Tr (AlE+i^E) -> G 32 . 

(A.ll) 

Following the convention of Ref . ^| for the twenty-seven 

127) = | 127,5) + i|27,i>. 

(A.12) 

We verify that G32 has the same Feynman rules as G b but with opposite sign in the 
case of vertices containing one Kaon and three pions. Thus, it follows that |27, h) 
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as denned in Eq. (|A. 6|) behaves as an octet (keeping in mind that the term Gb — G a 
behaves as an octet). Consequently the contribution of the |27) must be one-third 
of that of the octet for AI — 1/2 isospin amplitudes. This is only applicable for 
the imaginary part. In fact, our result for the imaginary parts in Table satisfies 
these identities, but it is not for the results produced in Ref. 1101 
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